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Abstract 

Using an angular momentum projected single particle basis, a pnQRPA approach is used to study 
the 2i'Pf3 properties of ten isotopes, exhibiting various quadrupole deformations. The mother and 
\ daughter nuclei exhibit different quadrupole deformations. Since the projected basis enables a 

unified description of deformed and spherical nuclei, situations where the nuclei involved in the 
^ . double beta decay process are both spherical, both deformed or one spherical and another deformed, 



can be treated through a sole formalism. Dependence of single /3 and strength distribution on 
atomic mass number and nuclear deformation is analyzed. For the double beta decay process, the 



■ Gamow- Teller transition amplitudes and half lives are calculated. Results are compared with the 

experimental data as well as with the predictions of other theoretical approaches. The agreement 
between the present results and experimental data is fairly good. 

PACS numbers: 23.40.Hc, 23.40.-s 
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I. INTRODUCTION 



The double beta decay of a (Z,N) nucleus, may take place through two distinct channels. 
For one channel, the final state consists of the residual nucleus (Z+2,N-2), two electrons and 
two antineutrinos, while for the other one the final state is lacking the antineutrinos. Sug- 
gestively, the two decay modes are called two neutrinos double beta (2z//5/5) and neutrinoless 
(Oz//5/5) double beta decay, respectively. The second mode is especially interesting since its 
discovery may provide a definite answer to the question whether neutrino is a Dirac or Ma- 
jorana particle. The 2i//5/3 process is interesting by its own but also due to the fact that 
may provide realistic nuclear matrix elements which might be further used to quantitatively 
evaluate the rate of the neutrinoless double beta decay. 

For such reasons many theoreticians focused their efforts in describing consistently the 
data for 2z//?/5 decay. The contributions over several decades have been reviewed by many 
authors. Instead of enumerating the main steps achieved toward improvingthe theoretical 



description we advise the reader to consult a few of the review works 

It is interesting to note that, most of the double beta emitters are deformed nuclei, the 
proposed formalisms use a single particle spherical basis. More than 10 years ago, one of 
us j9| proposed a formalism to describe the process of two neutrinos double beta decay in a 
projected spherical basis. It was for the first time that a proton-neutron quasiparticle random 
phase approximation (pnQRPA) for a two body interaction in the ph and pp channels with a 
deformed single particle basis was performed. Moreover, effects which are beyond pnQRPA, 
have been accounted for by means of a boson expansion procedure. Few years later the 
influence of nuclear deformation upon the contribution of the spin-flip conflgurations to the 
Gamow- Teller double beta transition amplitude, was studied 10|. In the meantime several 
papers have been devoted to the extension of the pnQRPA procedure to deformed nuclei, 
the applications being performed for studying the single beta decay properties as well as the 
double beta decay rates. Thus, pnQRPA approaches using as deformed single particle basis 
Nilsson or deformed Woods Saxon states have been formulated jlll. 12,[l^I3|- Also, a self- 
consistent deformed method was formulated where the single particle basis was obtained as 



eigenstates of a deformed mean fleld deflned throug 



1 a Hartree-Fock treatment of a density 
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dependent two body interaction of Skyrme type 

In a recent publication we continued the project opened in Ref.|2i|, by improving the 
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single particle basis. Indeed, in Ref.p] the single particle energies were depending linearly 
on a parameter which simulates the nuclear deformation. By contrast, in Ref. Q the 
core volume conservation constraint, ignored in the previous paper, determines a nonlin- 
ear deformation dependence for single particle energies. Of course, having different single 
particle energies, the pairing properties and the double beta matrix elements are expected 
to be modified. Another issue addressed in the previous paper was whether considering 
different deformations for the mother and daughter nuclei, modifies significantly the double 
beta transition amplitude (Mgt)- The answer to this question is positive since modifying 
the deformation for the daughter nucleus, the ground state correlations are modified and 
consequently the pnQRPA collapse point is changed. On the other hand the overlap matrix 
elements of the states describing the intermediate odd-odd nucleus, defined as excited states 
from the mother and daughter ground states respectively, are decreased. Therefore, consid- 
ering different nuclear deformations for mother and daughter nuclei quenches the Gamow- 
Teller (GT) double beta decay amplitude, which results in improving the agreement with 
the experimental data. 

The angular momentum projected spherical basis enables a unified description of spherical 
and deformed nuclei. Here we use this virtue of the single particle basis defined in Ref. jisl] 
and try to depict the specific features of the transitions between nuclei of similar or different 
nuclear shapes: a) spherical — > spherical (^^Ca ^^Ti), b) spherical ^deformed-prolate 
(^28Te ^^^Xe, ^=^°Te ^^^Xe), c) spherical -^deformed- oblate (^^^Xe ^ ^^^Ba, ^^^Xe 
^^^Ba), d) deformed — ^-spherical (^^°Pd ^^°Cd), e) deformed-prolate— >-deformed-prolate 
(962r^ 96mo), f) deformed-oblate^ deformed-oblate (^°°Mo ^ ^°°Ru, ^°^Ru -^^°^Pd, ^^^Cd 
— > ^^^Sn). It is worth mentioning that except for ^°^Ru, ^^'^Pd and ^^"^Xe for all other cases 
experimental data are available. 

The results of the present paper will be described according to the following plan. For the 
sake of a self-sustaining presentation, a brief review of the projected spherical single particle 
basis will be presented in Section II. Also, the basic equations necessary for calculating the 
GT double beta transition amplitude are given. In Section III, we discuss the results for ten 
double beta emitters: ^^Ca, ^^Zr, ^^^Mo, i°^Ru, "°Pd, "^Cd, ^^STe, ^^^Te, ^^^Xe, ^^sXe, for 
which the strength distribution for single (3~ emission, the Mqt and half lives values for the 
double beta decay process, are presented. Also, the strength distribution for the /5+ decay of 
the corresponding daughter nuclei is presented as function of the pnQRPA energy. A short 
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summary and concluding remarks are given in Section IV. 



II. PNQRPA TREATMENT OF THE GT (3(3 TRANSITION AMPLITUDE 
A. Projected single particle basis 

In Ref. jl^, one of us, (A.A.R.), introduced an angular momentum projected single 
particle basis which seems to be appropriate for the description of the single particle motion 
in a deformed mean field generated by the particle-core interaction. This single particle 
basis has been used to study the collective Ml states in deformed nuclei 171 as well as the 
rate of double beta process . Recently, a new version has been proposed where the 

deformation dependence of single particle energies is nonlinear and therefore more realistic 
jisl The new single particle basis has been used to study the double beta decay of 
deformed nuclei |l5]- In order to fix the necessary notations and moreover for the sake of a 
self-contained presentation, we describe briefiy the main ideas underlying the construction 
of the projected single particle basis. 

The single particle mean field is determined by a particle-core Hamiltonian: 

H = H,r^ + H,,,,-Muy E "v^Am- (2.1) 

A=0,2 -A</x<A 

where if^m denotes the spherical shell model Hamiltonian, while Hcore is a harmonic 
quadrupole boson (6^) Hamiltonian, associated to a phenomenological core. The interaction 
of the two subsystems is accounted for by the third term of the above equation, written in 
terms of the shape coordinates aoo, a2Af The quadrupole shape coordinates are related to 
the quadrupole boson operators by the canonical transformation: 

«2M = ;^(&L + (-r&2,-/.), (2-2) 

where k is an arbitrary C number. The monopole shape coordinate is to be determined from 
the volume conservation condition. 

Averaging H on a given eigenstate of Hgrnj denoted as usual by \nljm), one obtains a 
deformed quadrupole boson Hamiltonian whose eigenstate is an axially symmetric coherent 
state: 

^g = exp[d{bto-b2o)]\0)b, (2.3) 
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with \0)b standing for the vacuum state of the boson operators and d a real parameter 
which simulates the nuclear deformation. On the other hand, averaging H on \E'g one 
obtains a single particle mean field operator for the single particle motion, similar to the 
Nilsson Hamiltonian. Concluding, by averaging on a factor state of the particle core space 
the rotational symmetry is broken and the mean field mentioned above may generate, by 
diagonalization, a deformed basis for treating the many body interacting systems. However, 
this standard procedure is tedious since the final many body states should be projected over 
angular momentum. 

Our procedure defines first a spherical basis for the particle-core system, by projecting 
out the angular momentum from the deformed state 

"^Zj = \nljm)^,. (2.4) 

The upper index appearing in the l.h. side of the above equation suggests that the product 
function is associated to the particle-core system. The projected states are obtained, in the 
usual manner, by acting on these deformed states with the projection operator 

Pmk = / Di,^*{Q)R{Q)dQ. (2.5) 

We consider the subset of projected states : 

Kijid)=KijPLi[\nljI)^,]^K,Kfj{d). (2.6) 

which are orthonormalized and form a basis for the particle-core system. The main properties 
of these projected spherical states are: 

a) They are orthogonal with respect to I and M quantum numbers. 

b) Although the projected states are associated to the particle-core system, they can be 
used as a single particle basis. Indeed, when a matrix element of a particle like operator is 
calculated, the integration on the core collective coordinates is performed first, which results 
in obtaining a final factorized expression: one factor carries the dependence on deformation 
and one is a spherical shell model matrix element. Thus, the role of the core component 
is to induce a quadrupole deformation for the matrix elements of the operators acting on 
particle degrees of freedom. 

c) The connection between the nuclear deformation and the parameter d entering the 
definition of the coherent state ()2.3|) can be obtained by requiring that the strength of the 
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particle-core quadrupole-quadrupole interaction be identical to the Nilsson deformed term 
of the mean field. 

To the projected spherical states, one associates the 'deformed' single particle energies 
defined as average values of the particle-core Hamiltonian H' = H — Hcore- 



Since the core contribution to this average value, does not depend on the quantum num- 
bers of the single particle energy levels, it produces a constant shift for all energies. For this 
reason such a term is omitted in ()2.7|) . However, when the ground state energy variation 
against deformation is studied, this term must be included. 

In Ref.H H was shown that s.ngle p.tic.e ene.g,es. defined above, exh.bit a nonlinea. 
dependence on the deformation parameter d. Such a dependence is determined by the 
monopole-monopole interaction term after implementing the volume conservation constraint. 
Moreover, the deformation dependence of the new single particle energies is similar to that 
shown by the Nilsson .nodel Q. Therefore, the average values 4, may be viewed as 
approximate expressions for the single particle energies in deformed Nilsson orbits [20||. We 
may account for the deviations from the exact eigenvalues by considering, at a later stage 
when a specific treatment of the many body system is performed, the exact matrix elements 
of the two body interaction. 

Although the energy levels are similar to those of the Nilsson model, the quantum numbers 
in the two schemes are different. Indeed, here we generate from each j a multiplet of (2j + 1) 
states distinguished by the quantum number J, which plays the role of the Nilsson quantum 
number Q and runs from 1/2 to j. Moreover, the energies corresponding to the quantum 
numbers K and —K are equal to each other. On the other hand, for a given / there are 2J + 1 
degenerate sub-states while the Nilsson states are only double degenerate. As explained in 
Ref . , the redundancy problem can be solved by changing the normalization of the model 
functions: 



Due to this weighting factor the particle density function is providing the consistency result 
that the number of particles which can be distributed on the (2I-I-1) sub-states is at most 
2, which agrees with the Nilsson model. Here a stands for the set of shell model quantum 




(2.7) 




(2.8) 



M 
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numbers nlj. Due to this normalization, the states $q, used to calculate the matrix elements 
of a given operator should be multiplied with the weighting factor y^2/(2/ + 1). 

Concluding, the projected single particle basis is defined by Eq. ()2.(j|l . Although these 
states are associated to a particle-core system, they can be used as a single particle basis 
due to the properties mentioned above. The projected states might be thought of as eigen- 
states of an effective rotational invariant fermionic one-body Hamiltonian H^ff, with the 
corresponding energies given by Eq. ()2.7|) . 

HeffK^' = ei{d)^i''. (2.9) 

This definition should be supplemented by the request that the matrix elements of any 
operator between states and have, as we mentioned above, a factorized form, one 

factor carrying the d dependence, while the second one being a spherical shell model matrix 
element. Due to these features, these states can be used as single particle basis to treat many 
body Hamiltonians which involve one-body operators. This is the case of Hamiltonians with 
two body separable forces. As a matter of fact, such a type of Hamiltonian is used in the 
present paper. 

As shown in Ref. in the vibrational limit, d ^ 0, the projected spherical basis goes 
to the spherical shell model basis and e^;^ to the eigenvalues of Hsm- 

A fundamental result obtained in Ref.|3l for the product of two single particle states, 
which comprises a product of two core components, deserves to be mentioned. Therein we 
have proved that the matrix elements of a two body interaction corresponding to the present 
scheme are very close to the matrix elements corresponding to spherical states projected 
from a deformed product state with one factor as a product of two spherical single particle 
states, and a second factor consisting of a common collective core wave function. The small 
discrepancies of the two types of matrix elements could be washed out by using slightly 
different strengths for the two body interaction in the two methods. Due to this property 
the basis (2.6) might be used for studying any two-body interaction. 

B. The model Hamiltonian and its pnQRPA approach 

In the present work we aim at describing the Gamow- Teller two neutrino double beta 
decay processes with the property that mother and daughter nuclei may exhibit different 
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shapes. Indeed, in the chosen cases they might be both spherical, both deformed but with 
different deformations or one spherical and the other one of a deformed shape. The specific 
feature of the formalism used in the present work consists in treating all cases in an unified 
manner by using a sole single particle basis. 

The main ingredients of our formalism are as follows. The Fermi transitions contributing 
about 20% and the "forbidden" transitions are ignored, which is a reasonable approximation 
for the two neutrino double beta decay in medium and heavy nuclei. As usual, the 2vf3[5 
process is conceived as two successive single (3~ transitions. The first transition connects 
the ground state of the mother nucleus to a magnetic dipole state l"*" of the intermediate 
odd-odd nucleus which subsequently decays to the ground state of the daughter nucleus. 
The states mentioned above, involved in the 2v(3(3 process, are described in the framework 
of the pnQRPA formalism, by using the following many body Hamiltonian: 

+ 2xY.r,{pn)pt,{v'n'){-Y - 2xi E (2.10) 

The operator c^aihii'^raiM) creates (annihilates) a particle of type r (=p,n) in the state 
when acting on the vacuum state |0). In order to simplify the notations, hereafter the 
set of quantum numbers a{= nlj) will be omitted. The two body interaction consists of 
three terms, the pairing, the dipole-dipole particle-hole [ph) and the particle-particle [pp) 
interactions. The corresponding strengths are denoted by Grij = p,n),x,Xij respectively. 
All of them are separable interactions, with the factors defined by the following expressions: 



t _ ^ ^ J J 



M 



= E -Y{pIM\(Tf,\nl'M')^C^j^jMCnI'M', 
MM' ^ 



PiM) = E ^{pimcrM'M')^cl,,,c[j;^, (2.11) 

M,M' 

The remaining operators from Eq. ()2.1()|l can be obtained from the above defined operators, 
by hermitian conjugation. 

The one body term and the pairing interaction terms are treated first through the 
standard BCS formalism and consequently replaced by the quasiparticle one body term 
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J2tIM ^rO'tiM^TiM- In terms of quasiparticle creation (a^^^j^^) and annihilation (o^/m) oper- 
ators, related to the particle operators by means of the Bogoliubov-Valatin transformation, 
the two body interaction terms, involved in the model Hamiltonian, can be expressed just 
by replacing the operators ()2.1H) by their quasiparticle images which at their turn can be ex- 
pressed as linear combination of dipole two quasiparticle and quasiparticle density operators 
defined as: 

A'^ (nn) = filplnl „t „t 

^IfiKf"') / . rrin ^J. plpmp"'nlnmn^ 

mp,mn 

Bl^ipn) = ^ C'mp -m„ ^Ctp/pmp'^n/„m„( — )'^" "^^ = — [Op7p«^]lAc (2.12) 

mp,mn 

The quasiparticle Hamiltonian is further treated within the pnQRPA formalism, i.e. one 
determines the operator 

= T.[X{k)Al{k) - Y{k)A,_,{k){-y-% (2.13) 

k 

which satisfies the restrictions: 

[riM, = [H,p, rg = tuTl. (2.14) 

These operator equations yield a set of algebraic equations for the X (usually called 
forward going) and Y (named back-going) amplitudes: 

Z\\x(k)\' - lym'] = 1. (2.16) 

The analytical expressions for the pnQRPA matrices A and B are given in Ref.j3|. Since 
the pp interaction has an attractive character, for a critical value of xi the lowest root of 
the pnQRPA equations may become imaginary. Suppose that Xi is smaller than its critical 
value and therefore all RPA solutions (i.e. u) are real numbers and ordered as: 

< ^^2 < ••• < ^Ns- (2-17) 

Here Ng stands for the total number of the proton-neutron pair states whose angular mo- 
menta can couple to l"*" and moreover their quantum numbers n, I are the same. Hereafter 



the phonon amplitudes X and Y will be accompanied by a lower index "i" suggesting that 
they correspond to the energy cui. 

Since our single particle basis states depend on the deformation parameter d, so do 
the pnQRPA energies and amplitudes. The pnQRPA ground state (the vacuum state of the 
pnQRPA phonon operator) describes an even-even system which might be either the mother 
or the daughter nucleus. In the two cases the gauge and nuclear deformation properties are 
different which results in determining distinct pnQRPA phonon operators acting on different 
vacua describing the mother and daughter ground states, respectively. Therefore, one needs 
an additional index distinguishing the phonon operators of the mother and daughter nuclei. 
The single phonon states arc defined by the equations: 

Mj-^ki,\0)j, j^ij; 1,2,... N,. (2.18) 

Here the indices i and / stand for initial (mother) and final (daughter) nuclei, respectively. 
This equation defines two sets of non-orthogonal states, {|lfc/x)j} and {|lfcAi) /} describing the 
neighboring odd-odd nucleus. The states of the first set may be fed by a beta minus decay 
of the ground state of the mother nucleus while the states of the second set are populated 
with a beta plus transition operator from the ground state of the daughter nucleus. 

If the energy carried by leptons in the intermediate state is approximated by the sum of 
the rest energy of the emitted electron and half the Q- value of the double beta decay process 

AE^^Qpp + mec'', (2.19) 

the reciprocal value of the 2z//?/? half life can be factorized as: 

(T,7f )-^ = F\MGT{Ot ^ 0+)|^ (2.20) 

where F is an integral on the phase space, independent of the nuclear structure, while Mqt 
stands for the Gamow- Teller transition amplitude and has the expression : 

"^^"^^'S E, + AE + E,, ■ ^'-'^^ 

In the above equation, the denominator consists of three terms: a) AE, which was already 
defined, b) the average value of the k-th pnQRPA energy normalized to the particular value 
corresponding to k=l, i.e. 
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Ek = -j{,'^i,k + ^^f,k) - + (2.22) 

and c) the experimental energy for the lowest 1^ state. The indices carried by the transition 
operators indicate that they act in the space spanned by the pnQRPA states associated to 
the initial (i) or final (f) nucleus. The overlap matrix elements of the single phonon states in 
the mother and daughter nuclei respectively, are calculated within the pnQRPA approach 
and have the expressions: 

^{lk\lk')f = E W.pn)Xk'{f.pn) - Yk{i,pn)Yk{f,pn)] . (2.23) 

pn 

Throughout this paper, the Rose [2^ convention for the Wigner Eckart theorem is used. 

Before closing this section we would like to say a few words about what is specific to 
our formalism. Since our single particle states are projected spherical states, the pnQRPA 
formalism is fully identical to the one, usually employed for spherical nuclei. Since in the 
vibrational limit, {d ^ 0), our basis goes to the spherical shell model basis, one may say 
that the present formalism provides a unified description of spherical and deformed nuclei. 



III. NUMERICAL RESULTS 



The formalism described in the previous section, has been applied to the following ten 
isotopes: ^^Ca, ^^Zi, ^^^Mo, i°^Ru, ^^^Pd, "^cd, ^^^Te, ^^^Te, ^^^Xe, ^^^Xe. The spherical 
shell model parameters for these double beta emitters and the corresponding daughter nuclei 
are given by: 

hujo = AlA^/^, C = 2huoK, D = hujon, (3.1) 

with the strength parameters k and /i having the same (Z,N) dependence as in Ref. j^]. 

The angular momentum projected basis depends on two additional parameters. These are 
the deformation parameter d and the factor k entering the canonical transformation relating 
the quadrupole coordinate and boson operators. They are fixed in the same manner as in our 
previous publication [l^. Indeed, we require that the relative energy for the states 
and llfifl) be equal to that of Nilsson levels with Q = ^ and Q = ^ in the = 3 major 
shell. Moreover, adding to the mean field term defined before a QQ two body interaction 
we require that the lowest root for the charge conserving QRPA equation be equal to the 
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experimental energy of the lowest 2'^ state in the mother nucleus. Throughout this paper, 
the M-degenerate states are denoted by |n + 1 Ijl). 

The BCS calculation has been performed within a restricted single particle space. Due to 
the level crossing, the restriction of the single particle space for deformed nuclei is different 
from that for spherical nuclei. Indeed, in spherical nuclei Ikeda sum rule (ISR) is satisfied if 
two major shells plus the spin orbit partner of the intruder state are included in the single 
particle space. Suppose that the neutron open shell has N=3 with the intruder state |lg'9/2), 
in the standard spherical shell model picture. In the present formalism, including the spin- 
orbit partner state \lgl /2) means to consider the states ^q^^? with I ~ 7/2, 5/2, 3/2, 1/2. 
However, some of these states are higher in energy than states belonging to the |2d5/2/) 
multiplct. Due to such features appearing both in the upper part of the major open shell 
of neutrons and the bottom side of the proton major open shell we truncated the space 
considering an inert (Z,N) core and a number of states lying above the core states. The core 
and the number of outside states are chosen such that the non-occupation probabilities for 
the neglected bottom states as well as the occupation probabihties for the ignored upper 
states are smaller than 0.01. Of course, the single particle space for protons and neutrons 
are the same. Our calculations were performed with the core and number of states given in 
Table I. Once the single particle space is defined, the number of the dipole proton-neutron 
states can be calculated. Furthermore, the dimensions of the pnQRPA matrices for mother 
{Di) and daughter {D2) nuclei are readily obtained. These dimensions are also given in 
Table I. It is worth mentioning that using the single particle spaces given in Table I, Ikeda 
sum rule is satisfied for both the mother and daughter nuclei considered in the present paper. 



Nucleus 


48Ca 96zr iooMo 


lo^Ru 


llOpd 


ii6cd i28Te 




i34xe i36xe 


The (Z,N) core 


(0,0) (20,20) (26,26) (26,26) (26,26) (26,26) (44,44) (44,44) (44,44) (44,44) 


Number of states 


19 20 20 


22 


23 


27 22 


23 


21 23 


Di 


118 128 132 


140 


154 


166 142 


150 


138 154 


Dz 


115 128 132 


110 


151 


IGG 128 


132 


120 110 



TABLE I: The number of single particle proton states lying above the (Z,N) core is given. The 
single particle space for neutrons is identical to that for protons. Di and D2 are the dimensions of 
the pnQRPA matrix for mother and daughter nuclei, respectively 
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Note that despite the fact that single particle energies have a deformation dependence, 
we keep calling a major shell a set of states characterized by the same quantum number 
A^(= 2n + /) plus the states from the shell + 1 of maximum j. 

Single particle parameters d and k as well as the pairing strengths, fixed so that the mass 
difference of the neighboring even-even nuclei are reproduced, are listed in Table III. 

Now, let us turn our attention to the proton-neutron dipole interactions. In Ref[ll] it 
was suggested a simple A-dependence for these interaction strengths: 

X = ^MeV, xi = ^Mer. (3.2) 

We recall that this A dependence of the proton-neutron ph interaction strengths was obtained 
by fitting the position of the GT resonance for ^''Ca, ^°Zr and ^''^Pb. The pp interaction 
strength given above has been obtained by fitting the half lives for Z < 40 nuclei against 
the single l3~^ decay. A certain caution, however, is necessary when these formulae are used, 
since the A dependence is conditioned by t 



26|,|27|. For example, in Ref. 



le^mass region [2^ as well as by the single 
the GT resonance centroids in ^^^Te and 



particle space 

^^°Te, located at 13.7 and 14.1 MeV respectively, are reproduced with the x values equal to 
0.157 and 0.16 MeV respectively. These values for x are different from the predictions of 
Eq. ()3.2j) corresponding to A=128 and A=130, respectively. Moreover, as we see from Table 
IV,in the current paper the right position of these GT resonances are obtained by using 
X = 0.268 for both isotopes. It is noteworthy the fact that the daughter nuclei involved in 
a double beta process are stable against transitions. Therefore Xi is to be determined 
either using information about the half life of a l3~^ emitter lying close, in the nuclide chart, 
to the daughter nucleus under consideration or by fitting the data for a (p,n) reaction having 
the daughter as a residual nucleus. Hereafter, the ratio x/Xi is denoted, as usual, by gpp. 

The adopted procedure to fix the proton-neutron dipole interaction strengths is as follows. 
Whenever, in the intermediate odd-odd nucleus, the position of the GT resonance centroid is 
known, the ph interaction strength is fixed so that the above mentioned data is reproduced. 
As shown in Table II, for some of the isotopes considered in the present paper the log ft 
values associated to the (3^ /EC and j3~ transitions of the corresponding intermediate nuclei, 
are experimentally known. For these particular cases, x ^^id gpp are fixed by fitting the two 
mentioned experimental data. The log ft values were calculated by using the following 
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Mother 


Transition 


Intermediate 


Transition 


Daughter 


nucleus 


log ft 


nucleus 




nucleus 


iooMo 


13+ /EC 


iooTc 




iooRu 


Exp. 






4.66 




Th. 


4.61 




4.66 




lO^Ru 


[3+ /EC 






104pd 


Exp. 


4.32 




4.55 ^) 




Th. 


4.20 




4.62 




llOpd 


/3+/£;(7 


iiOAg 




iiOCd 


Exp. 


4.08 ^) 




4.66 ^) 




Th. 


3.86 




4.83 






/3+/£;c 








Exp. 


4.39t[5:|5 




4.662 




Th. 


4.05 




4.670 




i28Te 


13+ /EC 


128j 




i28Xe 


Exp. 


5.049 ^) 




6.061 ^) 




Th. 


4.930 




6.226 





TABLE II: The experimental and theoretical log ft values characterizing the (3^ /EC and (3 
processes of the intermediate nucleus ground state (l"*"). Experimental data are from: [4^ 

'=)Q/fl,^^Q. ^^B. ^^B> 



expression for ft: 

ft 6160 

Here |liM) denotes the first dipole phonon state in the intermediate odd-odd nucleus while 
|0) is the pnQRPA ground state. The low index "/" may take the value "z" and "/" 
depending whether the end state of the transition is characterizing the double beta mother 
or daughter nucleus. Therefore / = / is associated to single (3~ transition, while I = i to 
the (3^ / EC process. We chose qa = 1-0 in order to take account of the effect of distant 
states responsible for the "missing strength" in the giant GT resonance [Ij. For ''^Ca, we 
considered first x and Qpp (second row of Table IV) as given by Eq. ()3.2|l . To see the effect 
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of gpp on Mgt we repeated the calculations by keeping the same x as before but taking 
gpp = 0(third row of Table IV). It seems that fixing x as to reproduce the GT resonance 
centroid and taking Qpp = yields a better agreement with the experimental data. This 
situation is presented in the first row of Table IV. For ^^Zr, x was fixed by fitting the energy 
for the GT resonance centroid, while Qpp was taken as required by Eg 13.21 For ^^°Te we took 
the same x aiid Qpp as for ^^^Te. It is interesting to note that for this value of x the position 
of the GT resonance, at 14.1 Mev is nicely reproduced. As for the last two double beta 
emitters included in Table IV, there are available data neither for the GT resonance nor for 
the log ft values characterizing the single f3~ and j3^ / EC transitions of the corresponding 
intermediate odd-odd nuclei. For these isotopes we supposed for x and Qpp a similar linear 
1/A dependence as for ^'^°Te. 

The strength parameters x and Qpp determined in the manner described above are col- 
lected in Table III. They are also listed for each isotope in the first row of Table IV. These 
parameters yield double beta half-lives which are to be compared with the corresponding 
experimental data. The same parameters are used to calculate the single beta strength dis- 
tributions, shown in Figs 1-4. However, in order to have a fair comparison of the present 
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39|, in the second row of Table IV we give the 



results and those of Klapdor et al. 
results obtained with x and xi given by Eq. ()3.2|) . 

Once the parameters involved in the model Hamiltonian are fixed, the BCS and pnQRPA 
equations can be solved and the results be used in Eq. ()2.21|) to calculate the Mgt amplitude. 
Further, Eq. ()2.2()|) is used to calculate the half life of the 2i//5/? process. The phase factor F 
is not depending on the nuclear state structure and therefore was taken as in Refs. 
The values for F, used in this paper, correspond to qa = 1.254 (see the comments at the 
end of this section). Results for Mqt and T1/2 are given in Table IV. Therein one may find 
also the available experimental data. 
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Nucleus 


d 


k 


Gp [MeV] 


Gn [MeV] 


X [McV] 


gpp 


i\Qp0 + rUecA [MeV] 




0.3 


10.00 


0.65 


0.45 


0.180 


0.0 


2.646 


48Ti 


0.05 


2.00 


0.46 


0.36 


0.180 


0.0 




96Zr 


1.5 


10.20 


0.26 


0.26 


0.5 


0.112 


2.186 


96Mo 


1.2 


7.20 


0.3 


0.3 


0.5 


0.112 






-1.4 


10.00 


0.28 


0.26 


0.060 


1.600 


2.026 


iooRu 


-0.6 


3.6 


0.285 


0.220 


0.060 


1.600 






-1.55 


8.80 


0.26 


0.2 


0.150 


2.750 


1.161 


104p^ 


-1.35 


6.94 


0.26 


0.180 


0.150 


2.750 




llOpd 


-1.6 


6.00 


0.30 


0.32 


0.148 


2.450 


1.516 




-0.8 


3.06 


0.30 


0.18 


0.148 


2.450 






-1.8 


3.00 


0.20 


0.245 


0.238 


1.680 


1.916 




-1.2 


2.50 


0.18 


0.275 


0.238 


1.680 




i28Te 


0.5 


1.62 


0.27 


0.220 


0.268 


1.250 


0.946 




1.7 


6.50 


0.23 


0.220 


0.268 


1.250 




isoxe 


0.493 


1.88 


0.24 


0.210 


0.268 


1.300 


1.776 




1.4 


5.00 


0.24 


0.205 


0.268 


1.300 






-0.1 


1.95 


0.28 


0.300 


0.260 


1.261 


0.931 


i34Ba 


-0.468 


1.50 


0.24 


0.240 


0.260 


1.261 






-0.1 


1.80 


0.23 


0.29 


0.256 


1.243 


1.751 


i36Ba 


-0.698 


2.16 


0.19 


0.20 


0.256 


1.243 





TABLE III: The pairing and Gamow Teller ph interaction strengths are given in units of MeV. 
The ratio of the two dipole interaction ( particle-hole and particle-particle) strengths, denoted by 
gpp^ is also given. The list of the deformation parameter d and the factor k of the transformation 
(2.2) are also presented. The manner in which these parameters were fixed is explained in the text. 



16 



2i'P(3 decay 


[MeV] 


gpp 


IMgtI 
[MeV-1] 


T 1/2 [yr] 


present 


exp. 


Suhonen et al. 


Klapdor et al. 


48 Ca ^48xi 


0.180 
0.346 
0.346 


0.0 
0.112 
0.0 


0.043 
0.032 
0.036 


5.23 • 10^9 
9.27 • 10^9 
7.48 • 10^9 


(4.2±1.2) • 10^9 a) 




3.2 • 1019 1) 




0.500 
0.213 


0.112 
0.112 


0.113 
0.219 


1.66 • 10^9 
0.44 • 10^9 


(1.4+3-5) . ^ol9 a) 

V — U.5/ 


0.44 • 1020 2) 


5.2 • lOi'^ 


ioomo ^iooRu 


0.060 
0.207 


1.600 
0.112 


0.305 
0.212 


4.61 • 10^^ 
9.55 • 10^8 


(8.0 ±0.6) • 10^^ '^^ 
(0.115i°;J]3) . 1020 b) 
(0.033+0-02) . -^020 c,d) 


2.9 • lOis 3) 


1.8 • lOi^ 


104r^ ^104pd 


0.150 
0.201 


2.750 
0.112 


0.781 
0.343 


0.76 • 10^1 
3.95 • 10^1 






1.8 • 102i 
3.09 • 1022 6) 




0.148 
0.194 


2.45 
0.112 


0.263 
0.218 


15.85 • 10^9 
22.99 • 10^9 






5.0 • 10i9 
1.24 • 102i 6) 


ii*^Cd ^iisgj^ 


0.238 
0.187 


1.680 
0.112 


0.116 
0.069 


3.86 • 10^9 
10.96 • 10^9 


(3.2 ± 0.3) • 10^9 a) 


5.1 • 1019 5) 


8.3 • lOi^ 
3.75 • 1019 4) 


i28xe ^i28xe 


0.268 
0.174 


1.250 
0.112 


0.090 
0.127 


0.55 • 10^4 
0.28 • 10^4 


(7.2 ± 0.3) • 102"^ ^) 
(1.5 ±0.2) • 1024 ^) 


5.6 • 1023 5) 


1.2 • 1023 

5.7 • 1023 *) 


i30xe ^i30xe 


0.268 
0.172 


1.300 
0.112 


0.055 
0.091 


0.261 • 10^1 
0.097 • 10^1 


(1.5-2.8) • 1021 b) 
(2.7 ±0.1) • 102i ^) 
(0.75 ± 0.3) • 102i f) 


0.26 • 102i 5) 


1.9 • 10i9 
1.2 • 1020 *) 


134xe ^134Ba 


0.260 
0.169 


1.261 
0.112 


0.039 
0.040 


3.75 • 10^4 
3.49 • 10^4 






5.1 • 1022 

2.5 • 1023 *) 


136xe ^i36Ba 


0.256 
0.167 


1.243 
0.112 


0.039 
0.068 


5.102 • 10^0 
1.69 • 10^0 


> 8.1 • 1020 a) 


1.3 • 1020 5) 


6.0 • 1019 
3.3 • 1019 *) 



TABLE IV: The Gamow- Teller amplitude for the 2u(3(3 decay, in units of MeV~i, and the cor- 
responding half life {T1/2), in units of yr, are listed for ten ground to ground transitions. The 
experimental half lives for the transitions of ^^Ca(^)Ref. 
Ref. 



130 



Te Refs. 



2M 



b) Refs. 



231, 1) Ref. 122 



Ref. 



iooMo ( 

24(, b) RefM Ref.M, "'^ Ref.y|)), n^Cd (^) Ref.M), i28Te (^) Ref.|24|, ^) Ref. 3)> ) 



column the results reported in Refs.[l| 2), 

with the theoretical results from the last column reported in Refs 
(6)^. and 



i36Xe Refs.bj) are also given. In the second last 

n 

3J] 3) and _35] are given. Comparison is also possible 
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3a] (unmarked), 







The Darameters v and o„„ are also eiven. 



Before discussing the results presented in Table IV, we would like to discuss the strength 
distribution for single (3~ and transitions of mother and daughter nuclei, respectively. 
Thus, in Figs. 1-4 the strengths B'{GT)„ and B'{GT)^ for mother and daughter nuclei 
respectively , folded with a gaussian having the width equal to 1 MeV, are plotted as 
function of pnQRPA energies. These are equal to one third of the (3~ and (3^ strengths 
respectively, in the standard definition. Thus, the difference between the total strengths 
B'{GT)^ and B'(GT)j^, characterizing the mother nucleus, is to be compared with the sum 
rule (N-Z). The results of our calculations are to be compared with the available data for 
the GT giant resonance and single beta strengths given in 

Refs.0,0Q- 

At a glance 

one may see that while for Te and Xe isotopes most of the strength is concentrated in a 
narrow resonance, for the remaining nuclei the GT resonaces have a complex structure being 
spread over a large energy interval. Actually this feature is in agreement with experimental 



data showing in ^^^I and ^^°I a single peak at 13.7 and 14.1 MeV respectively 5l| while in 



^°°Tc and "^In |ij two peaks at (7.8,13.2) and (8.8, 14.30) MeV. The location of the above 
mentioned peaks are reasonable well reproduced. In the case of ^^''Mo, the first peak is 
higher than that centered at 12.3 MeV. Increasing the value of the repulsive ph interaction 
X, the ordering of the two peaks magnitudes is changed. 

The (3~ and (3^ strengths of '^'^Ca have been studied in Refs.js^ and 5^. Thus, the 
GT resonance has been populated in the reaction ^^Ca(p,n)'^'^Sc. This resonance is spread 
over an energy interval between 4.5 and 14.5 MeV. The result presented for B'{GT)- in 
Fig.l upper left panel, agrees with the experimental data. The total 3~ strength quenched 
with a factor of 0.6 55j, accounting for the polarization effects on the single beta transition 
operator, ignored in the present paper, is compared with the corresponding data in Table 
V. As shown therein, the agreement between the calculated strength and the corresponding 
data is reasonably good. The only known data for the total f3'^ GT strength is for ^^Ti: 

B{GT)+ = 1.42 ± 0.2 (3.4) 

Our calculations, corresponding to the first row of Table IV, predict for this strength the 
value 2.59. 

Comparing the f3~ strength distribution among 2qp states with those corresponding to 
pnQRPA states, one may conclude that the quasiparticle correlations accounted for by the 
pnQRPA approach, favors the displacement of the strength toward higher energies. This, 
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in fact, is due to the repulsive character of the ph interaction. As shown in Figs. 1-4, this 
effect is more pronounced for ^^Zr, ^^s.isorjg i34,i36^g_ 

From Table VI, it results that most of the f3~ strength is due to the transitions relating 
the proton and neutron g states. By contrast in ^^Ca, ^^Zr and ^^^Cd the single particle 
decays involving / states carry most of the transition strength. Also the GT resonance peak 
in ^^°Te is determined mainly by the transition in the 2d state. While in the lighter isotopes 
the transitions ul nl' where either / and /' are equal to | or | prevail, in Te and Xe 
isotopes the transitions z/| ti^ and z/| tc^ are dominant. 

Concerning the strength distribution shown in Figs. 1-4, right panels, the following 
features are to be noticed. The magnitude of this strength is much smaller than that of (3~ 
shown in the left panels. Moreover, the final states in the (3^ process are lying in the lower 
part of the spectrum, below 7.5 MeV. This suggests that the pp interaction may strongly 
influence the strength distribution among these states. The sensitivity of Z?"*" decay rate 
against the pp interaction was first noticed in Ref.Q]. Due to this feature the (3 (3 transition 
amplitude is also significantly affected by increasing Qpp. Since the ph and pp interactions 
are of different nature, one is repulsive and the other one attractive, one expects that the 
two interactions have opposite effects on the f3^ strength. When the pp interaction is large, 
comparing it with the ph interaction, the f3~^ strength in the quasiparticle picture is shifted 
toward the lower states. These are the cases of ^"^Pd, ^^''Cd, ^^^Sn. When both the ph and 
pp interactions are large, the strength of 2qp states are very much suppressed in the 
pnQRPA approach. Such situations are met for ^^^'^^°Xe, and ^^^'^^^Ba. 

By inspecting Table VIII, we conclude that the largest /J"*" strength is carried by the single 
particle proton- neutron transition in the shells Ig (for ^°°Ru,^°^Pd, ^^^Sn, ^^'^Xe, ^'^^Ba), 2d 
(for 96Mo, ii°Cd), If (for ^^Ti), 2f(for ^^'^Ba) and lh(in ^^oxe). Identifying the common 
2qp configurations carrying most efficiently the f3~ and f3~^ strengths from Table VI and 
VIII respectively, one may conclude which pnQRPA states excited from the mother and 
daughter ground state respectively, do maximally overlap each other and therefore bring 
large contribution to the double beta decay. In some of the depicted cases the state excited 
by the transition operator belongs to the GT resonance. An excellent example on this 
line is that of the transition '^^Ca^^^Ti, where the state at 10.326 excited from the ground 
state of ^^Ca, and that with energy of 5.957, excited from ^^Ti, have maximal overlap due 
to the 2qp state vr(3p|i), z/(3p||) which, in both /3~ and transitions, carries a large 
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strength. 

Let us now focus our attention on the GT double beta transition amphtude. This was 
calculated by means of Eq. (|2.21|) . where the energy shifts (j2.19p are those listed in Table 
III and the measured values for 1+ are collected in Table X. The states, energies and over- 
lap matrix elements involved in the Mgt expression, were calculated within the pnQRPA 
approach. The results corresponding to various sets of proton-neutron dipole interactions, 
fixed in the manner explained before, are listed in Table IV. Therein, the half lives of the 
2v[5f3 process are also given. The agreement with the available data is fairly good. 

Comparing our results with those of Klapdor et al. 0, 3| one may say that the half 
lives predicted by the present paper with the dipole interaction strengths given by Eq. (|3.2p 
are, without any exception, larger by a factor ranging from 2 (^°'^Mo) to 31 (^^Zr). Note 
that projecting out the gauge symmetry the results for Te isotopes are close to those given 
here for low values of Qpp. Comparing the results corresponding to x and Qpp fixed by fitting 
either the GT resonance centroid energy or the log ft value for (3^ / EC transition of the 
odd-odd nucleus, and the log ft value of the [3~ decay ending with the daughter nucleus, 
with those obtained with a renormalised pnQRPA equations and an adjusted Woods Saxon 
single particle mean field, we note that they are close to each other. 

Note that for ^^Ca the sets of (x, Qp-p)-, listed in the last two rows, provide half lives larger 
than experimental data, which suggests that the value of x niust be smaller than required 
by Eq. ()3.2|1 . Indeed, decreasing x to the value given in the first row of Table IV (=1.80), 
the agreement between the calculated GT resonance energy and the measured GT centroids 
could be improved. Moreover, the agreement with experimental data concerning T1/2 is also 
improved. By comparison one can see that the agreement quality obtained in the present 
paper is similar to that yielded by a full shell model calculation in Ref.js^. For this set 
of dipole interaction strengths, the quenched total strengths of (3~ transition of '^^Ca and 
transition of '^^Ti are equal to 15.65 and 2.59 respectively. The dominant peaks in the 
(3~ distribution correspond to the pnQRPA energies of 6.63 and 12.61 MeV. The carried 
strengths are 1.154 and 3.344 respectively. 

An interesting feature for the decay of ^^Ca was pointed out by the shell model studies 
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6^,1631. This refers to the cumulative effect brought by the low lying states in Sc, 
which actually yield the bulk contribution to the matrix element. The higher 1"*" states have 
a coherent destructive effect on the matrix element. It is worth investigating these aspects 
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within the present formahsm. Indeed, in Fig. 5 we plotted the transition amphtude Mqt as 
a function of the upper hmit of energies included in the defining equation (2.21). In other 
words, for a given E the energy Ef; defined by Eq. (2.22) and involved in the Mgt expression, 
is restricted by < E. We note that in five energy intervals this function is a monotonically 
increasing function of E, while in the following interval the transition amplitude decreases 
when states of higher energy are added. Also one notices a saturation effect, namely the 
contribution of states with energy larger than 16 MeV is very small. One may conclude 
that our results concerning the behavior of the double beta transition amplitude, are on a 
par with those of the shell model calculations. It is remarkable the fact that the maxima of 
Mgt{E) and (3^ strength are reached for similar energies. 
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Nucleus 


48Ca 96zr 


iooMo 


104Ru iiopd ii^cd 


i28Te 


i30Te i34xe 


i36Xe 


0.6E-B(GT)_ 


15.650 28.886 


30.040 


29.527 33.172 38.051 


43.340 


47.059 47.028 


50.703 






26.690 


32.700 


40.080 


45.900 





TABLE V: Total strengths for the Gamow- Teller (3 transition (first row) are compared with the 
available experimental data (second row) . Theoretical results are quenched with a factor of 0.6. 



Data for ^^^Mo and "'^Cd are from Ref. 



1W( 



52 1 while those for I28,i30'pg ^^.^ from Ref. 51]. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 




Transition 


Strength 


Transition 


Strength 


Transition 


Strength 




K3/ii) 


^7r(3/||) 


1.155 


z.(34i) 


^7r(34i) 


0.519 


K3/ii) 


^7r(34i) 
^7r(3/n) 


0.700 
3.344 


9621- 




-^(3/ii) 


0.602 


K3/ii) 


^7r(34|) 


1.040 


K3/ii 


-^(3/ii) 


5.842 






^7r(44|) 


0.467 


K3/|i) 

i/(3pii) 


-vr(3/4) 
^7r(3p||) 


2.473 

0.786 


K3/4) 


-7r(3/|i) 


2.212 






--(4^11) 


4.950 


K44i) 


^^(441) 


1.456 


K44I) 


^7r(44|) 


0.702 






-^(441) 


0.428 








K3/ii) 


^7r(3/|i) 


0.716 




1/(441) 


^7r(44i) 


1.005 
















K4d||) 


-444I) 


1.342 


K44i) 


-vr(44f) 


0.885 


K44I) 


- ^(441) 


0.563 




K44i) 


- ^(44i) 


2.722 


K44I) 


-7r(44|) 


0.400 


K3/|i) 


-^(3/I|) 


1.703 








1.321 








K3/M) 


^7r(3/|i) 


0.266 


llOp^j 


^(45ii) 


^7r(45|i) 


2.554 


K44i) 


^7r(44|) 


1.028 


-(44 i) 


^7r(44|) 


1.570 




K44i) 


- ^(441) 


1.080 


K44i) 


->7r(44i) 


1.171 


K3/|i) 


-7r(34i) 


0.294 




K44I) 


- ^(441) 


2.705 
















K44I) 


- vr(44|) 


1.668 


K44i) 


^7r(44|) 


0.769 


K3/4) 


-7r(3/|i) 


1.121 




j.(4d|i) 


^7r(44i) 


0.462 




-^(44i) 


0.927 


K3/ii) 


-^(3/H) 


2.647 














1.217 




- mil) 


1.971 


128r|.g 


K4ffi5) 


- ^(4ffii) 


0.446 


K44i) 


^7r(44|) 


0.467 


K44I) 


- vr(44|) 


12.483 




K45ii) 


- ^(441) 


0.380 


K44i) 


-^(441) 


0.198 


K44I) 


-7r(44|) 


7.316 


Te 


K45il) 


-^(441) 


0.462 


K5/i^i 


)-vr(54|) 


0.430 


K44I) 


-^^(441) 


2.543 




K4dii) 


-^(441) 


0.342 


-(441) 


^7r(44i) 


0.291 


K44i) 


-^(441) 


19.205 








0.713 


K44i) 


-^(441) 


0.659 


^(4d|i) 


^7r(44i) 


0.784 










i.(44i) 


^7r(44|) 


0.625 


K44i) 


- ^(441) 


21.050 


i36Xe 


K4ffii) 


- ^(441) 


0.710 


K44i) 


- vr(441) 


0.392 


K44I) 


-7r(44|) 


23.661 




i/(4d||) 


^7r(44|) 


0.321 




^7r(44i) 


0.709 









TABLE VI: The strengths carried by the pnQRPA states contributing to the first, second and 
third (if any) peaks from the upper panels of Figs. 1-4, are listed. On the left side of these numbers, 
the 2qp configurations closest in energy to the corresponding pnQRPA states, arc given. This is 
the dominant configuration of the chosen pn phonon state. The states ^^fj (see Eq. (2.6) are 
specified by the quantum numbers (Nljl) where N=2n-|-1. Also, the orbital angular momentum 
values 0,1,2, ..are mentioned by the letters s,p,d,.., respectively. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 




pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


pnQRPA energy 


Strength 




6.633 


1.155 


7.959 


0.519 


10.326 


0.700 












12.611 


3.344 




7.077 


0.602 


10.316 


1.040 


12.416 


5.842 




7.367 


0.467 


11.633 


2.473 


13.125 


2.212 








11.901 


0.786 








5.625 


0.428 


9.282 


1.456 


11.678 


0.702 




5.790 


4.950 






12.104 


0.716 




6.452 


1.005 












5.181 


1.342 


9.028 


0.885 


11.296 


0.563 




5.444 


2.722 


9.224 


0.400 


11.670 


1.703 




6.327 


1.321 






11.900 


0.266 


llOpd 


3.470 


2.554 


10.319 


1.028 


12.641 


1.570 




5.053 


1.080 


11.067 


1.171 


12.792 


0.294 




6.339 


2.705 












2.359 


1.668 


6.042 


0.769 


13.236 


1.121 




3.221 


0.462 


6.378 


0.927 


13.346 


2.647 








7.083 


1.247 


13.407 


1.971 


i28Te 


5.339 


0.445 


10.173 


0.467 


13.713 


12.483 




6.880 


0.380 


10.401 


0.198 


14.048 


7.316 


i30Te 


6.553 


0.462 


10.129 


0.430 


13.652 


2.543 




7.9G5 


0.312 


10.351 


0.291 


1-1.107 


19.205 




4.188 


0.713 


7.403 


0.659 


12.217 


0.784 








7.761 


0.625 


14.865 


21.050 




7.545 


0.710 


11.437 


0.392 


15.359 


23.661 




7.902 


0.321 


12.424 


0.709 







TABLE VII: The energies of the pnQRPA states which give the largest strength contributions to 
the peaks in Figs. 1-4, left panels. The carried strengths are also given. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 


Transition 


Strength 


Transition 


Strength 


Transition 


Strength 


48Ti 


vr(3pii)^K34i) 


0.153 


7r(3/ii)^z.(3/||) 


0.307 


7,(3/11)^^(3^5 1) 


0.562 


"'-'Mo 


\ 11' V 11' 


0.1-16 


^(lr/^§) — //(lr/§§) 


0.025 






iooRu 


7r(44|)^K45ii) 


0.405 


7r(44|)^K44i) 
7r(4(/||) ^ J^(45^) 

^^11' ^^11' 


0.509 
0.362 






104pd 


7r(4ff||) ^ 


0.441 


^(4^/|l)^K45|i) 
7r(44i) ^ i^(4d||) 

\ 11/ V 11' 


0.247 
0.110 


vr(45|i)^K45ii) 


0.023 




7r(4dii)^K445) 


0.273 


7r(45|i)^K45ii) 

7r(44l) ^ i^(4d||) 

V 11' V 11' 


0.110 

0.325 


7r(45ii)^K45ii) 


0.027 






0.391 


^(5/H)-K5/ii) 

;r(4oi|) ^i/(4o||) 

\'^22-' \<='22' 


0.117 
0.927 


7r(44i)-z.(4d||) 


0.019 


i28Xe 


7r(44i)^K44i) 
7r(5/i^^) ^ z/(5/i^|) 


0.011 
0.011 


7r(45ii)^K44i) 
7^(4597)^ /4 77) 


0.020 
0.041 








^(441)^^(4^33) 

7r(5/iiii) ^ z/(5/i^|) 

\ 11' \ 11' 


0.020 
0.028 


vr(45|i)^K45ii) 


0.026 








vr(5/H)^K5/ii) 

^(4^1 ) ^ viMm 


0.106 

0.086 










i36Ba 


^(4511)^^(4^73) 

^(441)^^(4^7 1) 


0.120 
0.124 











TABLE VIII: The same as in Table VI but for the right panels of Figs. 1-4. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 


pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


48Ti 


5.957 


0.153 


6.940 


0.307 


7.565 


0.562 




3.3G1 


O.llG 


1.221 


0.025 






iooRu 


2.099 


0.405 


3.617 

4.437 


0.509 

0.362 






104pd 


0.863 


0.441 


4.265 
5.300 


0.247 
0.110 


10.725 


0.023 




1.817 


0.273 


4.405 
4.531 


0.110 
0.325 


6.873 


0.027 




1.727 


0.391 


5.718 
6.378 


0.117 
0.927 


8.461 


0.019 




2.554 
3.246 


0.011 
0.011 


5.025 
5.479 


0.041 
0.020 








3.356 
3.446 


0.020 
0.028 


4.674 


0.026 








3.756 

3.803 


0.106 

0.08G 










i36Ba 


3.242 
3.534 


0.120 
0.124 











TABLE IX: The same as in Table VII, but for the right panels of Figs. 1-4. 
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Nucleus 


48Sc 


96Nb 


iootc 




llO^g 




128j 


130j 






Ei+ [keV] 


338 


1116 














58 


85 


177 


177 



TABLE X: The experimental energies for the first 1"*" states in the intermediate odd-odd nuclei are 
given in units of keV. Data are taken from Refs.j 
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4e 



47 



4S 
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501 1 ■ The states 



in ^^Sc and ^^Nb, at 338 and 1116 keV respectively, are not assigned with angular momentum and 
parity. Here we, ad hoc, suppose that they have the angular momentum equal to one and a positive 
parity. For ^^^Cs, there is no available data for energy levels. For this case we adopted the same 
excitation energy for the state l"*" as in ^^^Cs. Also for I28,i30j ^j^g energies for 1"*" are the same as 
in Ref.[51| 



Nucleus 


SSD 


present 


Mgt 




Mgt 


h/2 


ioomo 


0.211 


5.860-10 19 


0.305 


2.82-10 19 


104Ru 


0.616 


7.493-10 21 


0.781 


4.655-10 21 


llOpd 


0.551 


2.208-10 20 


0.263 


9.694-10 20 


i^^Cd 


0.421 


1.780-10 19 


0.116 


23.63-10 19 


i28Te 


0.032 


26.950-10 24 


0.090 


3.38-10 24 



TABLE XL The Mqt and ti/2 values obtained with the single state dominance {SSD) hypothesis. 
For an easy comparison we give also the values obtained within the present formalism. By contrast 
to the ti/2 values given in Table IV, here the half lives correspond to = ^■ 
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Before closing this Section we would like to say a few more words about the adopted 
procedure for fixing the dipole proton-neutron strength parameters. For the sake of a unitary 
treatment, the half lives of all double beta decaying nuclei were calculated by taking for qa 
the value 1.254. However, in five of the situations considered the single beta properties 
for the intermediate odd-odd nuclei are determined by supposing an effective value {gA=^-) 
for the axial-vector coupling strength, which might simulate the contribution of the higher 
energy states. Thus, although the nuclear matrix elements as well as the proton-neutron 
interaction strengths are similar for double and single transitions, we considered that the two 
sets of properties mentioned above are influenced by different parts of the proton-neutron 
QRPA excitation spectrum. In this context it is worth mentioning that some time ago Abad 
et al. advanced the single state dominance hypothesis (SSDH) which asserts that for 
double beta decay processes where the intermediate odd-odd nucleus has the state l"*" as 
ground state, most of the contribution to the double beta matrix element comes from the 
the first intermediate dipole state. If that hypothesis works, then the double beta process is 
dominated by two virtual and succesive single /3~ transitions, one from the ground state of 
the mother nucleus to the ground state of the intermediate odd-odd nucleus, while the other 
one from there to the ground state of the daughter nucleus. In the meantime the SSDH 



has been considered by many authors j6f 



68, 69,0, 71, 72, 73, Q 



As shown in Table 



X, four odd-odd isotopes have indeed the first 1+ as ground state. Moreover, there are also 
available data for (3~ and (3^ / EC transitions of ^^°I. Therefore, for all five isotopes from 
Table II we checked the SSDH validity by keeping in Eq. (2.21) only the first term from the 
sum and considering the states overlap equal to unity. Also the ti/2 values are calculated by 
considering gA = 1. The results are compared with those obtained by summing up over all 
pnQRPA states, otherwise keeping (7^ = 1 in order to have a fair comparison. Results are 
listed in Table XI. From there one may conclude that indeed our calculations confirm the 
SSDH for ^''^Mo, ^''^Ru and '^^'^Pd, but not for the remaining two double beta nuclei, ^^^Cd 
and ^^'^Te. Most likely for these cases the summation in the expression of Mgt should be 
extended from one to few states. 

It is worthwhile to address the question how stable are these results against changing 
the dimension of the single particle basis. We checked this feature with a positive result. 
To be more concrete let us describe the modifications obtained for ^^°Pd. For this isotope 
we increased the dimension Di from 23 to 27, otherwise keeping the same parameters for 
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single particle states as before. We changed the pairing strengths in order to preserve the 
pairing properties, i.e. to have the gap parameters unchanged. The new {Gp, Gn) for mother 
and daughter are (0.281,0.271)MeV and (0.2795, 0.1665)MeV, respectively. The pnQRPA 
matrix has the dimension D2 = 186. The ISR value devates from N — Z hj 3%. The proton- 
neutron interaction strengths have been changed to iXiQpp) = (0.13735, 2. 4)MeV in order 
to keep the log ft values for the single beta transitions of ^^°Ag close to the experimental 
data. The results for these observables are 4.84 for /5~ and 3.70 for (3^ /EC. The double 
beta transition amplitude and the T1/2 obtained under the new circumstances, are 0.2626 
and 15.881-10^^?/r. One notes that these values are very close to those listed in Table IV. 
One may conclude that the results are stable against enlarging the single particle space and 
moreover our choice for Di is motivated by the fact that ISR is satisfied. 



IV. CONCLUSIONS 

In the previous sections we completed the project started in Ref. [3| by studying the 
2u(5j3 decay of another ten even-even nuclei exhibiting various shapes. In the chosen cases 
the mother and daughter nuclei have the following shapes: a) both are spherical, b) both 
deformed-prolate, c) both deformed-oblate, d) one spherical and another deformed-prolate, 
e) both are near spherical but prolate, f) both are near spherical but oblate. The defor- 

nn 

mations obtained for the ten isotopes are similar to those of Ref . |^ . In some cases 
these are different from nuclear deformations reported in Ref.j^^. Indeed, for example in 



the present paper as well as in Ref.j57|, the quadrupole deformation for ^°°Mo is negative 
while in Ref. 0i] this is positive. Moreover, as shown in Ref. a negative deformation is 
reached in ^°^Mo. The reason for this discrepancy might be the fact that there the stationary 
points of the energy function are obtained in the space of quadrupole and hexadecapole de- 
formations while here only the quadrupole variable is considered. Moreover, here an angular 
momentum projected single particle basis is used. 

It is manifest the fact that an oblate to oblate single beta transition is involving single 
particle configurations which are different from those appearing in a prolate to prolate tran- 
sition. Indeed, suppose that a certain number of nucleons are distributed alternatively in 
a prolate and an oblate single particle levels and that in the first case the Fermi level for 
neutrons is characterized by a small quantum number /. In this case the (3~ strength for the 
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prolate to prolate transition is carried by single particle dipole transitions between states of 
low / as well as of large / but originating from the upper shell. By contrary, in the oblate 
to oblate transitions, the privileged transitions are those relating neutron and proton single 
particle states with large / and those of small I from the upper shell. Such cases can be 
easily identified in Tables VI and VIII. Due to the feature mentioned above the strength frag- 
mentation is expected to be more pronounced in the oblate to oblate transitions. Actually 
such a situation is met for i°°Mo, ^^''Ru, ^^°Pd and ^^^Cd. 

The structure of the peaks seen in Figs 1,2 in the (3~ strength distribution of the nuclei 
mentioned before is as follows. In ^°°Mo the peaks are determined by transitions inside the 
shells Ig (1st peak) 2d (2nd peak) and Ig and If (3rd peak). For ^°^Ru the first two peaks are 
determined mainly by the transitions from the shell Ig while the third one by the transition 
from the shell If. In ^^°Pd the following shells are involved in the transitions contributing 
most to the three peaks: Ig (1st peak), Ig, 2d (2nd peak) and Ig, If (3rd peak). For ^^^Cd 
only one shell contributes most to any of the three dominant peaks: Ig (1st and 2nd peaks). 
If (3rd peak). 

Note that for Te and Xe isotopes, the /5~ strength is mainly concentrated in one pnQRPA 
state. These nuclei are almost spherical (Te isotopes are soft prolate while Xe isotopes are 
soft oblate). Moreover, in the daughter nuclei the nuclear deformation has the same sign as 
in the corresponding mother nuclei. In ^^^Te and ^^°Te the dominant single particle state 
np transitions are ^{Ag^^) 7r(4(yf||) and u{Ad^^ 7r(4(i||), respectively. In Xe isotopes 
the np single particle transition i/(45f||) 7r(4(7||) , prevails. 

The transition amplitudes Mgt and half lives for the 2v(3(3 process were calculated within 
the pnQRPA approach by using a projected spherical basis. The agreement with the avail- 
able data is quite good. The adopted fitting procedure for the pn dipole interaction strengths 
yield large values for Qpp, in several cases. These values are not far from the critical value, 
where the pnQRPA breaks down. It is an open question whether for these transitions a 
good agreement with the data would be possible by keeping a small Qpp, but accounting for 
higher pnQRPA effects. Inspecting the Table IV, one can judge not only on the agreement 
of the present results with the experimental data but also on the comparison between pre- 
dictions of different theoretical approaches. Indeed, the agreement with experimental data 
is reasonably good. Although they are based on different formalisms as well as different 
single particle states the present results and those of Suhonen et a/.|l[ 34, 35 1 are not far 
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from each other. Comparing the results of the present paper, obtained with x s^nd Qpp given 
by Eq. ()3.2p and the corresponding predictions from Refs. jSy, |37|, |39|, one notices that they 
are quite different. 

It is worth mentioning that for ^°^Ru and ^^°Pd the pn dipole interactions are fully 
determined by fitting the data concerning the (3^ / EC and (3~ decay log ft values of the 
ground state (!'*') of the intermediate nuclei ^°^Rh and ^^°Ag, respectively. The predictions 
for the double beta emitter half lives are 0.76 ■ lO^^yr and 15.85 ■ lO^^yr. They are 40 and 



8 times smaller than the corresponding predictions of Ref . [37[ . Our prediction for the half- 
life of ^^^Xe, against double beta decay, is 3.75 ■ 10^^?/r which exceeds by a factor 15 the 
corresponding finding of Ref. 36. Issl]. 

The single state dominance is confirmed, by our formalism, to be valid for ^°°Mo, ^"'^Ru 
and ^i°Pd. 

Finally we may conclude that the projected spherical single particle basis provides a 
suitable framework for a unified description of the double beta properties of spherical and 
deformed nuclei. The results presented in the previous publication 3] and here constitute 
a good starting point for studying the higher pnQRPA contributions to the 2v(3j3 process 
as well as the transitions populating the daughter nuclei in an excited state. 
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FIG. 1: (Color on line) Single /3~ strength, for ^^Ca (upper-left panel), ^^Zr (middle-left panel), 
ioojyJq (bottom-left panel), and single (3^ strength for ^^Ti (upper-right panel), Mo (middle-right 
panel), ^'^'^Ru (bottom-right panel), folded with a Gaussian function having the width of 1 MeV, 
are plotted as a function of the energy within the BCS and pnQRPA approximation. The pnQRPA 
calculations correspond to the values of x ^-nd Qpp listed inside the graphs. 



35 



10 



o 

h 4 
2 




a 

h 4 

2 




o 

m 4 

2 




' 1 

_ j 

i 


' 1 ' 1 ' 1 

104-. ■ 

Ru" 


1 


-- 2qp 


1 
1 


pnQRPA 


1 

1 




=0.15 MeV - 


1 
1 
1 


1 

1 s 


=2.75 ~ 




1 

II 


p 


-A /] 


m 




r 




- 


[1 -- 2qp 




1 1 pnQRPA 




1 


=0.148 MeV- 






=2.45 ~ 


— / V 






~ j / 
/ 1 


1 1 1 


\ 

V 1 1 1 


_ 

- 


III 

, Cd - 




11 -2qp 




1 ; pnQRPA 




1 X= 
1 


D.238 MeV- 


- 


1 
1 
1 


g =1-68 


- 


1 
1 


«=pp _ 


1 

_ 1 

A 1 
/\ 1 


Im 


1 


- / \ / 
/ \-' I 

J 


, 1 , 


IV , 1 



5 10 15 20 
E [MeV] 



1.2 



0.9 



H 

S 0.6 



I ' I ' I ' I 



104 



— M 
1 1 
I I 
I I 
I I 
I I 
I I 
I I 
I I 
I 

I 



Pd 



2qp 

pnQRPA 
X=0.15MeV 




0.9 



H 

S 0.6 



110 




Cd 



2qp 

pnQRPA 
X=0.148 MeV 



g =2.45 



\ I 



116 



Sn 




— 2qp 

pnQRPA 

\^ X=0.238 MeV 



\ g =1.68 

1 "pp 



I . I 



5 10 15 20 
E [MeV] 



FIG. 2: (Color on line) Single f3~ strength, for ^'^^Ru (upper-left panel), ^^''Pd (middle-left panel), 
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FIG. 3: (Color on line) Single f3~ strength, for ^^^Te (upper-left panel), ^^^Te (middle-left panel), 
"^^^Xe (bottom-left panel), and single (3~^ strength for ^^^Xe (upper-right panel), ^^''Xe (middle-right 
panel), ^^^Ba (bottom-right panel), folded with a Gaussian function having the width of 1 MeV, 
are plotted as a function of the energy within the BCS and pnQRPA approximation. The pnQRPA 
calculations correspond to the values of x 9pp listed inside the graphs . 
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FIG. 4: (Color on line) Single f3~ strength, for ^^^Xe (left panel), and single f3~^ strength for ^^^Ba 
(right panel), folded with a Gaussian function having the width of 1 MeV, are plotted as a function 
of the energy within the BCS and pnQRPA approximation. The pnQRPA calculations correspond 
to the values of x and Qpp listed inside the graphs. 
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FIG. 5: Double beta transition amplitude Mqt, given by Eq.(2.21), is represented as function of 
energy E for ^^Ca. The summation over k, in EQ. (2.21), is restricted hy Ek < E where E^ is 
defined by Eq. (2.22). 



39 



